Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 January 2012 (MN I^TeX style file v2.2) 



Biparametric Adaptive Filter: detection of compact 
sources in complex microwave backgrounds 

M. Lopez-Caniego and P. Vielva 

Instituto de Fisica de Cantabria (CSIC-Universidad de Cantabria), Avda. de Los Castros S/N, 39005, Santander, Spain 



4 January 2012 



ABSTRACT 

In this article we consider the detection of compact sources in maps of the Cosmic 
Microwave Background radiation (CMB) following the philosophy behind the Mexi- 
can Hat Wavelet Family (MHWn) of linear filters. We present a new analytical filter, 
the Biparametric Adaptive Filter (BAF), that is able to adapt itself to the statistical 
properties of the background as well as to the profile of the compact sources, max- 
imizing the amplification and improving the detection process. We have tested the 
performance of this filter using realistic simulations of the microwave sky between 30 
and 857 GHz as observed by the Planck satellite, where complex backgrounds can be 
found. We demonstrate that doing a local analysis on fiat patches allows one to find a 
combination of the optimal scale of the filter R and the index of the filter g that will 
produce a global maximum in the amplification, enhancing the signal-to-noise ratio 
(SNR) of the detected sources in the filtered map and improving the total number of 
detections above a threshold. We conclude that the new filter is able to improve the 
overall performance of the MHW2, increasing the SNR of the detections and, there- 
fore, the number of detections above a 5<j threshold. The improvement of the new filter 
in terms of SNR is particularly important in the vicinity of the galactic plane and in 
the presence of strong galactic emission. Finally, we compare the sources detected by 
each method and find that the new filter is able to detect more new sources than the 
MHW2 at all frequencies and in clean regions of the sky. The BAF is also less affected 
by spurious detections, associated to compact structures in the vicinity of the galactic 
plane. 
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1 INTRODUCTION 

The emission of extragalactic point sources at microwave fre- 
quencies is known to be one of the most critical contaminants 
of the cosmic microwave background (CMB) anisotropics 
( de Zotti et al.|2005||Toffolatti et al.|2005| ). This foreground 
emission is a strong source of bias in the estimation of the 
CMB temperature and polarization angular power spectra. 
Even at the frequency range of 60 — 90 GHz (which covers 
an optimal window for observing the CMB), this effect is im- 
portant. In particular, point sources contaminate the power 
spect ra for multipoles t ^ 800 (Tegmark 1997} |Tucci et al 
2005), which is translated into a bias when determining the 
cosmological parameters from these quantities. In particular, 
the scalar spectral index n s and the optical depth r are two 
of the most biased parameter^] In addition, the foreground 
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emission of point sources also introduces a high level of non- 
Gaussianity in the CMB anisotropics (e.g., Argi ieso et al.| 
|2003| ). This non-Gaussian signal is a very important confu- 
sion noise when someone is probing non-standard models of 
structure formation, which actually predict a certain degree 
of non-Gaussianity on the CMB anisotropics, typically at a 
similar or lower level than the Gaussian deviation caused 
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by the point sources (e.g., Komatsu et al 
aT||2009 ). In addition to the role of the extragalactic point 



sources as CMB contaminants, characterising their proper- 
ties at the microwave frequencies is a very important field 
per se. For example, studying the number of objects per flux 
interval provides useful information to understand the his- 
tory of galaxy evolution (e.g., Gonzalez-Nuevo et al.||2008 



Tucci et al.|2011 ) 



A large number of works have been presented in the lit- 
erature to mitigate the impact of extragalactic point sources 
on the science that can be extracted from the analysis of 
the CMB anisotropics. Two global approaches are usually 
followed, being, in fact, complementary. On the one hand, 
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detection algorithms are applied to the CMB images for re- 
moving or masking the emission due to the brightest point 
sources. On the other hand, statistical modelling of the back- 
ground of the residual point sources is adopted. This mod- 
elling is followed in the determination of the cosmological 
parameters as well as when exploring the compatibility of 
the CMB signal with a given non-Gaussian scenario. This 
work lies in the former category: the detection of the bright- 
est point sources. The literature on this field is quite large, 
not only in relation to the development of specific tools for 
the CMB problem, but also to the adaptation of techniques 
originally developed in radio (e.g., [H5gbom| [l974) and op- 
tical (e.g., Bertin &; Arnouts||1996 ) astronomy. We refer to 



Herranz & Vielva (2010) for a complete and recent tutorial 



on the detection of compact sources in CMB maps. 

Some of the most popular tools to perform the de- 
tection and flux estimation of point sources are wavelets. 
Wavelets represent the simplest case of point source de- 
tection, since it is just based on a thresholding criterion, 
without any prior knowledge on the statistical properties or 
the background (as the Matched filter, e.g., Tegmark & de 



Oliveira-Costa 1998; Argiie so et al.|[2Q09] | Lopez- Caniego et 



al.| |2009 ~ or, even more, on the statistical properties of the 
background and the signal (as the Bayesian methods, e.g., 
Lopez-Caniego et al. 2005; Carvalho et al. 2009; Argiieso et 



aLl |2011| ). For CMB experiments with a PSF well defined 
by a Gaussian function, the Mexican Hat Wavelet (MHW, 
built as function proportional to the Laplacian of a Gaussian 



kernel) has proved to be a very good tool (e.g., Cay on et al.| 
2000). In addition, the wavelet scale R can be optimized to 
provide a larger number of detections, depending on the sta- 
tistical properties of the image ( |Vielva et al.|2001| ). Further 
improvement can be achieved by applying subsequently the 
Laplacian operator to the Gaussian kernel. This produces a 
series of wavelets known as the Mexican Hat Wavelet Family 
(MHWn, [Gonzalez-Nuevo et al.|20"Q6| where MHW1 is the 
standard MHW). The degree of derivation n is a reflection 
of the statistical properties of the background. In particular, 
low values of n define a filter that has most of the power at 
low frequencies (or equivalently, in the context of the CMB, 
at large scales). As the value of n increases, the power of the 
filter moves towards smaller scales. Therefore, if an image is 
dominated by large scale structures in the background, then 
one should expect a better point source detection efficiency 
for a higher value of n than for a smaller one. 

In |L6pez-Caniego et al.| ( [2006| , using simulations based 
on the Planck Reference Sky Model available at that time, 
it was shown that, on average, the MHW2 provides a larger 
number of point source detections than the MHW and the 
MF, and a better flux estimation (see |L6pez-Caniego et al7 



(2006) for further details of the simulations, and Tauber et 
al. (2010) for a description of the Planck mission). As a 



matter of fact, the MHW2 has been successfully applied to 
the WMAP maps, providing a larger number of detections 



Wright||2008n 2009) 



than the Matched filter used by WMAP ( |Hinshaw et al. 
2007| |Wright et al |2009| ) as well as other methods ( |Chen fe 



In this work we generalize the idea behind the MHWn, 
by defining a biparametric filter, ip g (R) where the parameter 
related to the background fluctuations g is allowed to vary 
in a continuous way. Therefore, the determination of the g 
and R parameters is done jointly, attending to the statistical 



properties of the background. The paper is organized as fol- 
lows. In Section 2 the biparametric filter is defined. We also 
discuss under which circumstances the new filter defaults 
to other known filters. The performance of the method is 
illustrated in Section 3 by analysing in detail some simu- 
lated microwave images. A statistical analysis on how the 
biparametric filter behaves as a function of the frequency 
and the galactic latitude is presented. Conclusions are given 
in Section 4. 



2 METHODOLOGY 

As mentioned above, considering the statistical properties 
of the background in the vicinity of a compact object is a 
key step in the design of a filter to be used for the detec- 
tion of point sources. This issue can be as simple as looking 
for the optimal scale that provides the maximum amplifica- 
tion for wavelets or as sophisticated as providing a complete 
likelihood and prior functions, for instance, when applying 
Bayesian approaches. In all-sky CMB experiments such as 
WMAP ( |Bennett et al.|2003| ) or Planck ( |Tauber et al.|2010 ) 
one can deal with very different backgrounds. For instance, 
the properties of the background have a typical dependence 
on the Galactic latitude: the major emission corresponds 
to the Galactic components (synchrotron, free-free, dust) 
near the Galactic plane, whereas the CMB and the back- 
ground of extragalactic sources are, typically, more impor- 
tant at intermediate and high latitudes. In addition, depend- 
ing on the observational frequency, radio emissions (galactic 
and extragalactic) are more important in the lower part of 
the microwave frequency range (i.e., from 10 to 100 GHz), 
whereas thermal dust and the cosmic infrared background 
are the major contributors to highest frequencies (i.e., from 
300 to 1000 GHz). However, at intermediate frequencies the 
CMB is, overall, the most important background. Each one 
of these emissions has its own specificities and, therefore, it 
is suboptimal to use the same filter at all frequencies and 
positions in the sky. 

In this paper we explore the detection of objects in dif- 
ferent types of backgrounds commonly found in CMB exper- 
iments. These objects are, in general, point-like sources that 
have been convolved with the beam point spread function 
(PSF) of the instrument used for the observation. They can 
be described as a signal s(x) = Ar(x), where A is the intrin- 
sic flux density of the object and t(x) is the PSF of the beam 
and x is a unit vector in the sky. Although the profile r{x) 
can be described by any function, the formalism adopted 
in this paper assumes that it is given by an isotropic PSF, 
i.e., r(x) = t(|x|) = r(x). We adopt this simplification be- 
cause for most of the CMB experiments an effective isotropic 
window function is usually defined, even when the PSF has 
a certain degree of anisotropy. In any case, of course, the 
filter proposed in this paper can be easily generalized for 
non-circular beam profiles. 



2.1 The Biparametric Adaptive filter 

Following the idea behind the MHW Family, in this paper 
we propose a filter that has two free parameters that will 
allow us to better incorporate into the filter the statistical 
properties of the background. As for the case of wavelets, 
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one of the parameters of the filter is the scale R. The fil- 
ter scale is associated with a compression/expansion of the 
typical scale of the PSF and provides the size of the filtering 
kernel. In addition, we incorporate the index of the filter g 
that can be seen as a generalization of the role played by 
the order of the Laplaci an operator used in the defin ition of 
the MHW Family (e.g., |Gonzalez-Nuevo etal . 20060 It is 
related to the filter location in Fourier space or, conversely, 
with the filter oscillations in the real space. As it will be 
discussed in the next subsection, this parameter is somehow 
associated with the shape of the angular power spectrum of 
the background. The filter definition is such that it is com- 
pensated (e.g. Cay on et al.|2 000), and, therefore, for the case 



adopted in this paper of 2D Euclidean images, it behaves as 
a wavelet. The only case in which there is no compensation 
and, therefore, the filter is not a wavelet, is for the par- 
ticular case in which the background is defined by a white 
noise field. In this particular situation, the filter defaults to 
a Gaussian kernel with g = (see next subsection). 

Let us denote the BAF (defined in terms of the two 
parameters R and g) by ^(x; R, g,b). It denotes the value 
of the filter centred at 6, in the position x. Following the 
philosophy behind the wavelets, this filter can be seen as 
the scaling version of a translated mother filter: 



V(x;R,g,b) = j^iPg 



R 



(1) 



This mother filter tp g defines the BAF and, in Fourier space, 
is given by: 



1>g (qR) 



(qR) 9 T(qR). 



(2) 



In particular, and following the convention in |Sanz, Herranz^] 
|fe Martinez- Gonzalez (2001), convolving such a filter with 
a 2D image f(x) we obtain the filter coefficients uj g (R,b) 
defined as: 



dx f(x)V(x;R,g, 



(3) 



A very common situation in astronomy is to have a PSF de- 



fined by a Gaussian function r(x) — (l/27r<7b ) e 



where Gh is the Gaussian beam dispersion. For this particu- 
lar case, the BAF is given by: 



(qR) 



(qR) 3 



(4) 



and the filter coefficients at the position of a source with a 
profile Ior(qR), where Iq is the amplitude of the source, is 
given by: 



9 + 2 

(l + « 2 )^ 



(5) 



where z = R/(Jb- As it was introduced by Vielva et al. ( 2001 ), 
the filter parameters are determined by imposing a maxi- 
mum amplification A of the point source amplitude in the 
filter coefficients map: 



2 Notice that, whereas the order n of the Laplacian operator is a 
natural number, the index g is defined as a real number ^ 



A. 



0Jg(R)/<J U 



(6) 



where a is the dispersion of the image in the real space, Io 
is the amplitude of the source and is the dispersion of 
the filter coefficients at the scale R and index g, that can 
be defined as: 



(7) 



al(R,g) = 27T / dq qP (q) $(Rq), 
Jo 

with P (q) being the angular power spectrum of the analysed 
image. 



2.2 The BAF and other filtering kernels 

The BAF introduced in equation [2] defaults, under certain 
conditions, to other well known filter kernels extensively 
used in literature. In this sense, the BAF can be consid- 
ered as a generalization of these filters. The most obvious 



generaliz ation is for the MHWn (e.g., |Gonzalez-Nuevo et 



|al.||2006[ ip n ). As previously mentioned, when the profile 
of the point sources is described by a Gaussian function, 
the BAF is given by equation [4] Both filters are related by: 
i\) g — (2 n /tv) ijj n , where g = 2n and n = 1, 2, 3, — 

Another interesting situation occurs when the back- 
ground can be described by an angular power spectrum fol- 
lowing a power law P(q) = Cq _1 , for q > 0, and C being a 
normalization constant. In this case, it is easy to show, from 
equation [7] that the variance of the filter coefficients is given 
by: 



al(R,g) = C 



7-2 "p 



( 2g + 2- 1 



,7-2 



(8) 



[r(W 

From equation [6] one can demonstrate that the maximum 
amplification occurs when R = Oh and g = 7, i.e., the BAF 
defaults to the Matched filter. 



3 THE SIMULATIONS 

In order to test the performance of the new filter we will 
use realistic simulations of the microwave sky at 30, 44, 70, 
100, 143, 217, 353, 545 and 857 GHz. These simulations 
have been generated with the pre-launch Planck Sky Model 
(Lea ch et al.|2 008; Delabrouille et al. 2012), a software pack- 
age developed within the Planck Collaboration that allows 
one to make simulations at the microwave frequencies of 
the CMB, Galactic diffuse emissions and compact sources. 
The Galactic interstellar emission adopted in this paper is 
described by a three component model of the interstellar 
medium comprising of free- free, synchrotron and dust emis- 
sions plus added small scale fluctuations to reproduce the 
non- Gaussian nature of the interstellar emission. Free- free 



emission is based on the model of Dickinson et al. (2003) 



assuming an electronic temperature of 7000 K, where the 
spatial structure of the emission is estimated using a Ha 
template corrected for dust extinction. Synchrotron emis- 
sion is based on an extrapolation of the 408 MHz map of 



Haslam et al. ( 1982 ) from which an estimate of the free- 
free emission was removed. The thermal emission from in- 
terstellar dust is estimated using model 7 of Finkbein er et| 
al. (1999). Point sources are modelled with two main cate- 



gories: radio and infrared. Simulated radio sources are based 



© 0000 RAS, MNRAS 000, 000-000 



4 M. Lopez- Caniego and P. Vielva 




+1.500E+04 



Figure 1. From top to bottom, a combination of noise, CMB, 
galactic diffuse emission and compact source emission maps at 
30, 143 and 857 GHz simulated with the pre-launch Planck Sky 
Model. 



on th e NVSS or SUMSS ( |Condon et aT|l998| |Mauch et al. 
2003t and GB6 o r PMN catalogues fl Griffith et al.|1995||Gre- 
gory et al.|1996 ). Measured fluxes at 1 and/or 4.85 GHz are 
extrapolated to Planck frequencies assuming a distribution 
in flat and steep populations. Infrared sources are based on 
the IRAS catalogue (|Beichman et al. 1988), and modelled 



as dusty galaxies. A detailed description of each component 
can be found in |Leach et al. (2008). These simulations use 
the HEALPix pixelization scheme ( Gorski et al.||2005| with 
NSIDE=1024 at 30, 44 and 70 GHz and NSIDE=2048 for 
the rest. In Figure [I] we show three of these simulations, at 
30, 143 and 857 GHz. 



3.1 Patch analysis 

As it was mentioned above, the purpose of this work is to 
study the performance of the new filter in different back- 
grounds. For this reason, and to better illustrate the prob- 
lem, we have studied three different regions of the sky at 
each of the following three frequencies, 30, 143 and 857 GHz. 
For each region, we have projected flat patches of 7.3 x 7.3 



square degrees. At 30 GHz each patch has 128 x 128 pixels 
and a pixel size of 3.43 arcminutes. At 143 and 857 GHz, 
each patch has 256 x 256 pixels and a pixel size of 1.71 
arcminutes. Each region has been selected by visual inspec- 
tion to have increasing background complexities. In general, 
the statistical properties of the background in those patches 
with lower Galactic latitudes are more complex than those 
in the higher latitudes. The coordinates of the centres of the 
patches can be found in Table [I] In Figure [2] one can see the 
nine selected regions. 



In order to determine the filter characteristics that best 
adapts to the properties of each patch we look for the opti- 
mal scale R and index g of the filter that maximize the am- 
plification of the sources present in the patch. An increase in 
the amplification will allow one to detect more sources above 
a certain threshold. In Table ^ one can see the values of g 
obtained for each patch and each frequency. Note that the 
optimal values of g and R vary from one region of the sky to 
another. For the 30 and 143 GHz cases, these variations are 
small even for different Galactic latitudes, particularly for 
the 143 GHz case, but at 857 GHz, where the Galactic dust 
and the far-infrared background dominate, the variations are 
large. This simple analysis already shows the importance of 
optimizing not only the scale but also the index of the filter. 

To further illustrate this exercise, in Figure [3] we show, 
for the selected regions, the amplification as a function of 
the scale R and index g of the filter. Since the behaviour of 
the amplification as a function of the optimal scale R has 

Vielva et al.] 



been already studied in previous works (e.g. 
2001 Lopez-Caniego et al.|2 006), we will concentrate on the 
properties of the index g. In this respect, in Figure ^] we 
show the results of the same analysis but representing the 
amplification as function of the index g, conditioned on the 
optimal scale R. As we mentioned above, for each region 
there is a combination of R and g that produces a single 
maximum in the amplification. In particular, if we look to 
the upper left panel of Figure ^] corresponding to the region 
1 at 30 GHz, the maximum in the amplification is reached 
for g — 6.66. If we use, instead, the standard MHW at the 
optimal scale , g — 2, the amplification is only 82% of that 
of the maximum. If we compare the MHW2 at the optimal 
scale , g — 4, the amplification is ~ 97%. If we look at one 
of the regions for the 143 GHz case, for example the middle 
central panel of Figure [2] the maximum appears at g — 6.5. 
In this case we can see a similar behaviour for the MHW 
at the optimal scale as before, the amplification is ~ 83% 
of that of the maximum and for the MHW2 at the optimal 
scale the amplification is ~ 97%. More interesting is the case 
of 857 GHz, where the maximum in g changes significantly 
from one region of the sky to another. In this frequency, the 
maximum of the lower left panel of Figure [4] is at g — 1.83 
and the MHW and MHW2 at the optimal scale produce 
an amplification of ~ 95% and ~ 100%, respectively. Then, 
the maximum in the lower central panel is at g — 4.33, 
and the MHW and MHW2 at the optimal scale produce an 
amplification of ~ 94% and ~ 100%, respectively. Finally, 
the maximum in the lower right panel is found at g — 6.66, 
and the MHW and MHW2 at the optimal scale produce an 
amplification of ~ 82% and ~ 97%, respectively. 
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Figure 2. Selected regions with increasing background complexities extracted from the 30, 143 and 857 GHz simulated maps. At 30 
GHz these patches have 128 x 128 pixels in size, where the pixel size is 3.43 arcminutes. At 143 and 857 GHz these patches have 256 x 256 
pixels in size, where the pixel size is 1.71 arcminutes. The upper panels correspond to the 30 GHz regions, the middle panels to the 143 
GHz regions and the lower panels to the 857 GHz regions. The left, center and right panels correspond to regions 1, 2 and 3 respectively. 
The coordinates of the center of these three regions can be found in Table ^ 



Region 


GLAT 


GLON 


#30 


#143 


#857 


1 


+60 


70 


6.66 


4.33 


1.83 


2 


-36 


290 


6.66 


6.50 


4.83 


3 


+01 


123 


5.66 


6.83 


6.66 



Table 1. This table shows the galactic coordinates of the centres 
of the three regions that we have selected visually in the simu- 
lated maps at 30, 143 and 857 GHz in order to have increasing 
background complexities. We also show the optimal value of the 
index of the filter g that we have found in each case. 



3.2 All-sky analysis 

In the previous section we have studied the performance of 
the BAF optimizing the two parameters, the scale of the fil- 
ter R and the index of the filter g, that produce a maximum 
in the amplification of the filtered map as compared with the 
unfiltered one. This study was done looking at three visually 
selected regions in the sky for each of the three considered 



frequencies. The results for this analysis with only nine re- 
gions already show that choosing the right index g and scale 
of the filter R can increase the performance of the filter sig- 
nificantly. In this section we want to improve the statistics 
of the analysis by increasing the percentage of sky and fre- 
quency coverage of the study. For this purpose, we use nine 
full-sky simulations between 30 and 857 GHz, dividing each 
of the nine simulations into 1344 flat patches that effectively 
overlap to cover the 100% of the maps. Then we apply our 
maximization techniques to obtain the pair of optimal val- 
ues R and g that define the filter, taking into account the 
local statistics of the background in each particular patch. 
In addition, we make considerations in terms of galactic lat- 
itudes since we know that the complexity of the background 
increases when one gets closer to the galactic plane. In Fig- 
ure [5] one can see the results of this analysis with respect to 
the index g for the nine frequency bands. In this figure we 
show the average filter index and its dispersion per Galactic 
latitude bin, where the bins have been chosen to cover the 
same area in the sky. 

For the 30 GHz case, the mean of the average filter in- 
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Figure 3. Point source amplification as a function of the index g vs. the scale z = R/cr^ of the filter for a selection of regions. The 
upper, middle and lower panels correspond to the 30, 143 and 857 GHz maps respectively. These panels correspond, from left to right, 
to the regions 1, 2 and 3. Every contour line represents l/25th of the maximum amplification. 



dex is rsj 7, except for the galactic region between [-20,20] 
where the synchrotron radiation dominates and the index 
g decreases a bit (see upper left panel of Figure [5}. Then, 
the behaviour of the filter index between 44 and 143 GHz 
is fairly flat, although it varies between 4-6. It is important 
to note that the behaviour of g starts to change again at 
217 GHz and becomes more and more steep when we get 
closer to 857 GHz, extending each time to higher galactic 
latitudes. This is a clear correlation between the contribu- 
tion of the dust emission and the Far-Infrared background 
to the background signal in the patches. It reflects how the 
filter is adapting to background complexity. This effect can 
be clearly seen in the lower panels of Figure [5] 

In order to do a qualitative comparison between the 
MHW2 and the BAF, we crossmatch the catalogs of objects 
detected with each technique at each frequency and use the 
common objects to build several sets of figures. 

First, in Figure [6] we represent the relative difference 
between the noise level estimated for one method and the 
other as a function of galactic latitude and frequency. This 
allows us to see the overall behaviour of the noise estimated 
with the new filter compared with the MHW2. In the up- 
per panels one can see that the estimation of the noise in 



the patches filtered with the BAF is a few percent smaller 
than that of the MHW2. When we increase the frequency 
one starts to see a change in the behaviour, as expected, 
and between 143 and 857 GHz the estimation of the noise 
of the BAF is up to a 25% smaller than that of the MHW2, 
but only in the vicinity of the galactic plane. Note that even 
though the BAF, by definition, includes the MHW2 and its 
performance should always be equal or better than that of 
the MHW2, in the 545 GHz panel of Figure [6] we see that 
for one of the bins, the relative difference is negative. This is 
not a problem of the filter but a rare artefact of the imple- 
mentation of the algorithm that occurs due to the fact that 
we divide the sky into overlapping patches and, sometimes, 
the same source is detected with the same technique in two 
adjacent patches but with different SNR. Since we keep the 
ones with highest SNR, it could happen that a source is de- 
tected with a higher SNR with the BAF in one patch and 
with the MHW2 in the adjacent one. In this case we cannot 
guarantee that the noise estimation of the BAF is equal or 
higher than the MHW2 because they are in fact looking at 
slightly different regions of the sky. 

Second, with respect to the flux density estimation of 
the sources obtained with the BAF or the MHW2 tech- 
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Figure 4. Point source amplification as a function of the as a function of the index g, conditioned on the optimal scale R. The upper, 
middle and lower panels correspond to the 30, 143 and 857 GHz maps respectively. These panels correspond, from left to right, to the 
regions 1, 2 and 3. Note that the amplification has been normalized to the maximal amplification in each case. 



niques, in Figure [7] one can see that between 30 and 217 
GHz the fluxes obtained with both methods follow a clear 
one-to-one line. At 143 GHz one starts to see an increase 
in the dispersion and a small bias that is more obvious at 
higher frequencies. These differences mitigate when all the 
objects that lie within a ±3 degrees cut in latitude are ex- 
cluded. These plots can be seen in Figure [8] where most of 
the scatter and bias has now disappeared. 

To further investigate this correlation between biased 
fluxes and extended galactic regions at high frequencies, we 
have plotted the flux of the MHW2 vs. that of the BAF for 
353, 545 and 857 GHz using an additional galactic cut of 
±15 degrees. In Figure [5] we show the fluxes without apply- 
ing any galactic cut (upper panels), applying a ±3 degree 
cut (middle panels) and applying a ±15 degree cut (lower 
panels). Again, the scatter and bias found in the upper pan- 
els decreases significantly when increasing the galactic cut, 
as we expected. 

This result is telling us that both niters produce unbi- 
ased flux density estimations outside the galactic plane and 
is highlighting a problem in its vicinity. We suspect that 
the problem may be related to any or both of the follow- 
ing effects. First, the fact that the extremely bright and 



spatially variable background signal is contributing to the 
recovered flux density of the point sources. In other words, 
the niters are not able to fully remove the contribution from 
the background. Second, the niters are not able to recover 
unbiased flux density estimations of the extended sources 
present in the galactic plane, because they have not been 
designed for this purpose and their response to these ob- 
jects can be different. In order to test these ideas, we have 
performed additional simulations injecting 200 point sources 
with the same flux density and spatially distributed along 
the galactic plane of the 545 GHz simulated map. Then we 
have attempted to estimate their flux densities using the 
BAF and MHW2 techniques and obtain that we are able to 
recover them with differences smaller than 1%, i.e., much 
lower than the ones in Figure [9] Therefore, we are confident 
that the reason for the bias seen at 217 GHz and above is 
due to the uneven performance of both niters when dealing 
with extended objects, a situation for which they have not 
been designed. 

Moreover, to illustrate the performance of the new filter 
we will do a simple comparison with the MHW2 in terms 
of the number of detections above SNR > 5 obtained by 
each method when applied to the nine simulated maps. In 
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Figure 5. The mean filter index g per the galactic latitude bin and its dispersion are shown. The galactic latitude bins have been chosen 
to cover an equal area on the sky. In blue the mean index g across the sky is given. 



the left and right panels of Figures |10[ |11| and [12] we show 
the positions in the sky of the sources detected by the BAF 
(left panels) and those detected by the BAF or the MHW2 
only (right panels). In these figures one can see that the 
number of objects detected in the vicinity of the galactic 
plane, as well as other complex regions such as the LMC, 
Orion, etc., is large and experience tell us that a fraction 
of those detections are not true sources but bright compact 
galactic structure that looks like sources. Starting at 217 
GHz, one can notice that the number of detections increases 
rapidly due to the change in the dominant population of 
sources, from non-thermal radio sources to thermal infra- 
red sources. In addition, one can see that the BAF seems to 
detect more new sources all over the sky as opposed to the 
MHW2 that concentrates most of its new detections in the 
vicinity of the galactic plane, many of which tend to be spu- 
rious detections caused by filaments and extended structures 
in the galaxy. This is an important result because it shows 
the potential of the BAF, a filter that better adapts itself 
to the local properties of the background, removing part of 
the noise and large scale emission more effectively. This will 
imply a reduction in the number of detections in complex 
regions (that are likely to be spurious) while increasing the 



number of real detections in cleaner regions of the sky due 
to the improved estimation of the noise. 

In Table |3.2| we show the number of detections with 
a SNR > 5 for both techniques when applied to the same 
regions of the sky, in one case covering the whole sky and 
in the other case applying a large galactic cut of ±30 de- 
grees. One can see that the BAF detects, in general, more 
sources, specially between 30 to 143 GHz. In some chan- 
nels the improvement is almost inexistent, for example at 
70 GHz. Here, due to the properties of the background, the 
optimal value of g is essentially 4, which corresponds to the 
MHW2. On the contrary, in the upper frequency channels 
(353, 545 and 857 GHz) the number of SNR > 5 sources 
produced by the MHW2 is a bit larger than those of the 
BAF. Note that these results apply when considering the 
detection above SNR > 5 in the whole sky, including the 
galactic plane where we know that the MHW2 is detect- 
ing more sources likely to be galactic emission rather than 
extragalactic compact sources. If one looks at the second 
case, where a galactic cut has been applied, the BAF de- 
tects more sources than the MHW2 in all cases except for 
one, 217 GHz, where the difference is very small. Note that 
these numbers are in agreement to what it was mentioned 
above. If we inspect the right panels of Figures [To] [TT] and 
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Figure 6. The average relative improvement (of BAF with respect to MHW2) in cr^ is given. The BAF is able to improve the estimation 
of the noise up to 25% in particularly complex regions, as in those areas dominated by galactic emission. 



[12] one can see that the MHW2 tends to detect more objects 
in complex regions than the new filter, many of which are 
most likely not true sources. This is telling us that the new 
filter is removing the background more efficiently which im- 
plies not only an improved estimation of the SNR but also 
a decrease in the number of spurious detections. 

In addition, we have compared the catalogs of detec- 
tions obtained with the BAF and the MHW2, searching 
for those objects detected only by one method and not by 
the other. We find that the BAF provides a higher num- 
ber of unique detections in the maps from 30 to 217 GHz, 
whereas the MHW2 detects more objects at the highest fre- 



quencies (see Table 3.2 ). A deeper comparison is made when 
distinguishing between unique detections in highly contami- 
nated regions (e.g., in the vicinity of the Galactic plane) and 
cleaner areas in the sky. Obviously, as mentioned above, the 
former will correspond (at high probability) with spurious 
detections associated to extended objects, whereas the lat- 
ter will correspond most probably with true point sources. 
In order to identify the contaminated regions, we generate 
two types of masks. The first one includes the 15% of the 
brightest pixels in each map. The second one, more conser- 
vative than the other, includes the 25% of brightest pixels in 
each map. We use the less conservative 15% masks between 



Table 2. This table shows the number of detections above 
SNR > 5 that we have found applying the MHW2 and the BAF 
techniques at the nine simulated maps. In addition we show the 
number of detections above SNR > 5 and galactic cut of ±30 de- 
grees to have an idea of the effect of the galactic emissions when 
doing this kind of simple comparisons. 
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Figure 7. Comparison of the estimated flux density of the common sources detected by both methods for channels 30 to 217 GHz. 





Figure 8. As in the previous figure, but applying a ±3 degrees three galactic cut to show that the small bias found in the 143 and 217 
GHz panels rapidly disappears when we exclude from the figure sources very close to the galactic plane 



30 and 217 GHz and the most conservative one between 
353 and 857 GHz, where the galactic dust emission extends 
to higher latitudes. In practice, these masks identify very 
well the galactic emission and a few complex regions across 
the sky (e.g., Magellanic clouds, Orion, Ophiucos). In Table 
|3.2| one can also see that BAF detects more unique objects 
both inside and outside the mask between 30 and 217 GHz, 
whereas in the most contaminated channels at 353 GHz and 
above, it also detects more unique objects outside the mask 



while the MHW2 detects up to three times more unique ob- 
jects inside the mask, many of which are most likely galactic 
extended emission rather than extragalactic point sources. 

Finally, in order to give an idea of the spatial distribu- 
tion in the sky of the values of the filter index g, the opti- 
mal scale R and the noise estimation of the filtered maps 
that we have obtained, we have constructed a set of figures 
assigning at each HEALPix pixel (NSIDE=8) the average 
value of g, R and of the detected sources that fall into it 
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Figure 9. Comparison of the estimated flux density of the common sources detected by both filters for channels 353, 545 and 857 (left, 
center and right panels) and applying three galactic cuts (none, ±3 degrees and ±15 degrees) to show that the bias found in the figures 
rapidly disappears when we exclude from the figure sources very close to the galactic plane 



(see Figures 13 14 and 15). It is interesting to note that, as 
expected, the values of the optimal index g decrease a bit 
in the galactic plane at 30 GHz and significantly increase in 
and around the galactic plane at 217 GHz and above, while 
they remain fairly homogeneous across the sky between 70 
and 143 GHz. A similar behaviour can be seen in the op- 
timal scale R showing that the scale of the filter increases 
with increasing complexity of the background. In addition, 
one can see that the estimated noise in the filtered maps a u 
follows the galaxy and a few other complex regions in the 
sky very well. 



4 CONCLUSIONS 

In this work we have developed and studied the performance 
of a new filter that maximizes the amplification of the com- 
pact sources embedded in a complex background using two 
free parameters, the scale R and the index g of the filter. 
This new filter is called biparametric adaptive filter (BAF). 
To study the capabilities of this new filter we have used sim- 
ulations of the microwave sky at the frequencies of Planck 
between 30 and 857 GHz. 



In a first detailed analysis of three frequencies (30, 143 
and 857 GHz) we have visually selected three regions of in- 
terest with increasing background complexity and obtained 
projected patches of 7.3 x 7.3 square degrees. We have ap- 
plied the BAF to these regions and have demonstrated that 
we can always find a combination of R and g that maximizes 
the amplification of the sources in the filtered map with re- 
spect to the original map. In addition, we have found that 
it is important to optimize not only the scale of the filter 
but also the index g that adapts to scaling properties of the 
background, finding values of g in the range [1.8-12] for the 
nine discrete patches that we have studied. 

In a second test, we have explored the performance of 
the BAF when doing a full-sky analysis. We divide the sky 
into 1344 overlapping patches, 7.3 x 7.3 square degrees each, 
that effectively cover the 100% of the sky and apply our new 
filtering technique to each one of them, detecting sources 
above SNR > 5. The results not only confirm what we 
found in the preliminary study of nine interesting regions 
but also show how the index g changes in a smooth and 
coherent way when we move from low galactic latitudes to 
high ones in the presence of strong galactic emission such as 
the synchrotron radiation at 30 GHz or the dust emission 
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Figure 10. Position of the detected sources at 30 (upper panels), 44 (medium panels) and 70 GHz (lower panels) for the two considered 
techniques, MHW2 and BAF. In the left panels we show the position in the sky of all the objects detected above SNR > 5. In the right 
panels we show only those objects SNR > 5 that were detected by one method and not by the other, and viceversa. 



and the Far-Infrared background starting at 217 GHz and 
above. The most extreme cases are 353, 545 and 857 GHz, 
where the index g changes from very low values close to 
g — 2 for the highest galactic latitudes up to g — 12 for the 
regions very close to the galactic plane. Even more, in these 
channels that we could consider cleaner from galactic emis- 
sion, in the sense that the CMB is the dominant component 
of the background, one can see that the behaviour of the in- 
dex is fairly flat at all galactic latitudes and takes values, on 



average, between g — 4 and 6. In particular, at 70 GHz, the 
index is always very close to g = 4, the index of the MHW2. 
In addition we have qualitatively compared the performance 
of the filter with that of the MHW2, in terms of the number 
of detections above SNR > 5 in two cases, in the whole sky 
and above a galactic cut of ±30 degrees. We find that in the 
first case, the BAF detects more sources than the MHW2 
below 217 GHz and the MHW2 detects more source at 217 
GHz and above, although from the inspection of the position 
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Figure 11. Position of the detected sources between 100 and 217 GHz for the two considered techniques, MHW2 and BAF. In the left 
panels we show the position in the sky of all the objects detected above SNR > 5. In the right panels we show only those objects 
SNR > 5 that were detected by one method and not by the other, and vice versa. 



of the sources in the sky one can see that most of the new 
detections at high frequency are very close to the galactic 
plane and are likely to be spurious detections due to bright 
compact emission from the galaxy. In the second case, were 
the galactic cut has been applied, we find the BAF detects 
more sources at all frequencies. Moreover, we have looked 
at the number of unique detections obtained by one method 
and by the other, inside and outside a galactic mask, and 
concluded that, first, the BAF detects more unique objects 



at all bands inside and outside the mask up to 217 GHz, 
and, second, that at 353 GHz and above the BAF detects 
less objects than the MHW2 inside the mask, where most of 
the detections are spurious, as mentioned above. 

We have demonstrated that a tool to detect compact 
sources like the MHW2 is a very good compromise for the 
kind of backgrounds that one can find in microwave exper- 
iments, both at low and high galactic latitudes, at all the 
Planck frequencies, but its performance can be improved if 
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Figure 12. Position of the detected sources between 353 and 857 GHz for the two considered techniques, MHW2 and BAF. In the left 
panels we show the position in the sky of all the objects detected above SNR > 5. In the right panels we show only those objects 
SNR > 5 that were detected by one method and not by the other, and viceversa. 



we use the BAF, a filter that explores the combination of 
the scale R and index g that best adapt to the profile of the 
sources and to the local properties of the background max- 
imizing the amplification and the SNR of the detections. 
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Figure 13. Spatial distribution across the sky of the values of the filter index g found during the analysis of nine simulations between 
30 and 857 GHz. The upper panels show the 30, 44 and 70 GHz cases, the middle panels show the 100, 143 and 217 GHz cases and the 
lower panels show the 353, 545 ad 857 cases. 




Figure 14. Spatial distribution across the sky of the values of the optimal scale of the filter z = R/(Jb found during the analysis of nine 
simulations between 30 and 857 GHz. The upper panels show the 30, 44 and 70 GHz cases, the middle panels show the 100, 143 and 217 
GHz cases and the lower panels show the 353, 545 ad 857 cases. 
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Figure 15. Spatial distribution across the sky of the estimated noise in the filtered maps (7^ found during the analysis of nine simulations 
between 30 and 857 GHz. The upper panels show the 30, 44 and 70 GHz cases, the middle panels show the 100, 143 and 217 GHz cases 
and the lower panels show the 353, 545 ad 857 cases. 
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Table 3. In this table we present the total number of detections 
above SNR > 5 obtained by one method and not by the other, as 
well as the number of detections inside and outside a a mask that 
we have defined for each frequency and that includes the 15% of 
the brightest pixels in the maps between 30 and 217 GHz and a 
more conservative 25% of the brightest pixels between 353 and 
857 GHz. We use these masks to define what can be considered 
as a complex region and what is not. 
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